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ABSTRACT 


This  thesis  presents  a  Computation  Fluid  Dynamics  (CFD)  analysis  of  the 
aerodynamic  perfonnance  of  circular  and  elliptical  leading  edges  of  compressor  blades, 
with  a  range  of  leading  edge  droop  angles.  Specifically,  simulations  were  conducted, 
with  a  free  stream  Mach  number  of  0.65  to  quantify  the  change  in  pressure  distributions 
and  boundary  layer  momentum  thickness  in  the  leading-edge  region  for  a  range  of 
incidences  on  a  flat  plate  with  various  leading-edge  ellipticity  ratios,  ranging  from  unity 
(circular)  to  5.5.  In  addition,  the  impact  of  drooping  the  leading  edge  was  analyzed  over 
a  range  of  incidence  angles  from  zero  to  13  degrees.  Pressure  distributions  indicate  that 
elliptical  leading  edges  can  eliminate  separation  bubbles  at  zero  incidence.  The  results 
indicated  that  the  minimum  loss  occurred  at  an  ellipticity  ratio  of  about  3.5  and  elliptical 
leading  edges  with  a  droop  slightly  greater  than  the  average  incidence  can  significantly 
decreased  aerodynamic  losses  over  a  wider  range  of  incidences. 
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I.  INTRODUCTION 


A.  OVERVIEW  AND  BACKGROUND 

The  United  States  Navy  depends  heavily  on  the  use  of  gas  powered  turbine 
engines  for  power  and  propulsion  for  most  of  it  marine  and  aviation  platforms.  However, 
with  the  increasing  cost  of  fossil  fuels,  and  the  ever  present  threat  of  tightening  defense 
budgets,  there  is  an  increasing  focus  on  making  these  power  plants  more  efficient.  To 
date,  thousands  of  studies  have  been  conducted  with  the  aim  of  increasing  the  overall 
efficiency  of  gas  turbines,  and,  more  specifically  to  increase  the  efficiency  of  the 
turbine’s  compressor.  Methods  of  reducing  losses  have  been  analyzed  through  the  intense 
study  of  the  aerodynamics  of  the  compressor’s  blades.  In  addition  to  the  military, 
industries  in  the  civilian  sector  stand  to  benefit  from  an  increase  in  overall  gas  turbine 
efficiency. 

One  way  to  influence  the  efficiency  of  the  compressor  is  to  modify  the  geometry 
of  its  blades  and  hence  their  aerodynamic  characteristics.  In  the  1950’s  and  1960’s,  the 
NACA  65  series  blade  became  popular  for  use  in  axial  compressors  as  an  innovative  and 
more  efficient  blade  type.  Figure  1  is  a  cross-section  of  a  NACA  65-213  airfoil. 


Suction  Side 


Flow  Direction  Pressure  Side 

Figure  1.  NACA  65-213  Airfoil  Section 

The  geometry  of  any  airfoil  can  be  broken  down  into  three  main  parts.  They  are 
the  leading  edge  shown  on  the  left  end  of  Figure  1,  the  trailing  edge  shown  on  the  right 
end,  and  the  body  of  the  airfoil  connecting  the  leading  and  trailing  edges.  The  lifting  side, 
shown  as  the  top  in  Figure  1,  is  referred  to  as  the  suction  side,  and  the  bottom  side  is 
known  as  the  pressure  side. 
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The  most  important  parts  of  compressor  blade  geometry  are  the  leading  edge  and 
the  trailing  edge.  The  camber  and  thickness  distributions  of  the  body  are  important,  but  as 
long  as  it  smoothly  connects  the  leading  and  trailing  edges,  its  geometry  has  secondary 
impacts  on  the  flow  when  compared  to  the  leading-edge  direction  and  shape  and  trailing 
edge  direction. 

Today,  many  gas  turbines  use  a  circular  shaped  leading  edge,  as  can  be  found  on 
the  LM2500  gas  turbine,  for  example.  The  circular  leading  edge  shape  is  popular  because 
it  is  easy  to  manufacture,  and  maintains  its  shape  relatively  well  in  the  hostile 
environments  encountered.  Figure  2  is  a  close  up  view  of  the  nose  of  the  NACA  0012 
airfoil,  which  has  a  circular  shaped  leading  edge. 


Figure  2.  Circular  Leading  Edge  of  a  NACA  0012  Airfoil  (see  Ref.  23) 

Air  flowing  around  the  leading  edge  of  the  blade  plays  an  important  role  in 
determining  the  aerodynamic  characteristics  of  the  flow  over  the  rest  of  the  blade.  As  the 
air  approaches  the  leading  edge  of  the  blade,  it  stagnates,  and  from  there  accelerates 
around  the  blade,  going  over  either  the  suction  side  or  the  pressure  side.  The  accelerated 
flow  then  may  or  may  not  undergo  a  rapid  deceleration,  depending  on  the  current 
incidence  and  other  flow  factors.  If  the  geometry  of  the  leading  edge  and  the  flow 
conditions  result  in  a  large  enough  deceleration,  the  flow  adjacent  to  the  surface  of  the 
blade  may  experience  such  an  adverse  pressure  gradient  that  it  actually  travels 
backwards.  Figure  3  depicts  the  velocity  vectors  of  an  attached  boundary  layer  flow, 
whereas  the  diagram  in  Figure  4  shows  flow  that  has  separated. 
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Figure  3.  Attached  Boundary  Layer  Velocity  Profile 


Figure  4. 


Separated  Boundary  Layer  Velocity  Profile  75 


If  separation  occurs,  the  flow  can  either  reattach  along  the  surface  of  the  blade,  or 
remain  unattached.  A  region  where  the  flow  separates  and  then  quickly  reattaches  is 
called  a  separation  bubble.  The  length  of  these  separation  bubbles  is  usually  quite  small; 
on  average  less  than  five  percent  of  the  total  chord  length.  Separation  bubbles  most 
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commonly  occur  very  close  to  the  leading  edge  where  the  flow  undergoes  such  dramatic 
acceleration. 

The  goal  of  many  compressor  blade  designers  has  been  to  create  a  geometry  that 
will  allow  the  airfoil  to  operate  without  the  generation  of  a  separation  bubble.  In  cases 
where  the  fonnation  of  a  separation  bubble  is  imminent,  mitigating  their  effects  can 
dramatically  improve  the  overall  performance  of  the  blade.  However,  comparatively  little 
work  has  been  done  with  a  focus  on  controlling  the  flow  through  manipulation  of  the 
leading-edge  geometry. 

One  way  to  decrease  the  chances  of  forming  a  separation  bubble  on  the 
compressor  blade  is  to  redesign  the  shape  and  direction  of  the  leading  edge.  When  the 
flow  goes  over  a  circular  leading  edge,  there  is  very  rapid  acceleration.  However,  shapes 
that  allow  for  the  flow  to  more  gradually  accelerate  and  decelerate  have  been  shown  to 
reduce  the  likelihood  of  an  adverse  pressure  gradient  causing  a  separation  bubble  from 
occurring.  The  elongated  nature  of  an  ellipse  has  been  shown  to  better  control  the 
acceleration  and  deceleration  of  the  flow  around  the  leading  edge. 

An  ellipse  is  generally  defined  as  the  distance  from  loci  lying  on  the  semimajor 
axis.  The  relative  distance  between  the  loci  is  described  using  an  equation  relating  its 
semimajor  and  semiminor  axes,  and  is  referred  to  as  eccentricity.  The  ratio  of  the 
semimajor  axis  to  semiminor  axis  is  referred  to  as  the  ellipticity  ratio.  The  greater  this 
ratio,  the  more  elongated  the  ellipse.  Figure  5  below  depicts  a  comparison  of  an  elliptic 
leading  edge  to  a  circular  leading  edge. 
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Figure  5.  Comparison  of  Elliptic  and  Circular  Leading  Edges  Attached  to  a  Flat  Plate 

In  Figure  5,  an  elliptic  leading  edge  is  compared  to  a  circular  leading  edge,  both 
attached  to  the  same  airfoil  body.  The  circular  leading  edge’s  semimajor  and  semiminor 
axes  are  the  same  length,  resulting  in  an  ellipticity  ratio  of  unity.  The  elliptic  leading  edge 
has  a  greater  semimajor  axis  is  noticeably  less  blunt  when  compared  to  the  circular 
leading  edge.  This  geometric  alteration  of  the  forward  section  of  the  blade  allows  the 
flow  to  accelerate  and  decelerate  around  the  leading  edge  much  more  gradually. 

In  addition  to  the  leading  edge  shape,  the  leading  edge  direction  plays  an  equally 
important  role  on  impacting  the  flow  and  aerodynamic  performance  of  an  airfoil.  One 
way  researches  have  observed  the  impact  of  the  leading-edge  direction  is  by  drooping  the 
forward  section  of  the  airfoil.  The  acceleration  of  the  airflow  around  the  leading  edge 
becomes  greater  with  increasing  incidence,  and  drooping  the  leading  edge  of  the  airfoil 
allows  the  nose  to  point  more  into  the  oncoming  airflow  at  these  higher  incidences.  In 
this  way,  the  airfoil  can  perfonn  as  an  undrooped  leading  edge  would  at  lower  incidences, 
and  allow  the  blade  to  perform  better  at  higher  incidences.  Figure  6  shows  a  comparison 
of  a  drooped  elliptic  leading  edge  to  that  of  a  basic  elliptic  leading  edge. 
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Figure  6.  Comparison  of  Straight  Elliptic  Leading  Edge  and  Elliptic  Leading  Edge  with  a 

Droop  Angle  8 

B.  LITERATURE  SURVEY 

Several  researches  have  conducted  studies  of  the  impact  of  leading-edge  shape 
and  direction  on  the  aerodynamic  characteristics  of  flow  over  various  airfoils.  Their 
experimental  and  computational  work  supports  the  theory  that  leading-edge  ellipticity 
and  droop  can  significantly  impact  pressure  distributions  and  aerodynamic  losses. 

1.  Previous  Studies 

The  experiments  of  Walraevens  and  Cumpsty  [1]  compared  the  aerodynamic 
characteristics  of  a  circular  leading  edge  attached  to  a  flat  plate  to  those  of  a  flat  plate 
with  a  leading-edge  ellipticity  ratio  of  1.89.  In  their  studies,  testing  conditions  were  such 
that  the  flow  remained  nearly  incompressible.  They  witnessed  that  the  leading-edge 
geometry  directly  affects  the  aerodynamic  performance  of  the  airfoil  through  its  impact 
on  formation  of  a  separation  bubble.  Figure  7  shows  a  comparison  of  the  pressure 
distributions  over  the  two  different  airfoils. 
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Figure  7.  Walraevens-Cumpsty  Pressure  Coefficient  (-Cp)  as  a  Function  of  Distance  Along 

the  Plate  From  the  Geometric  Leading  Edge 

Their  investigation  revealed  that  for  the  circular  leading  edge,  the  pressure  spike 
reached  much  higher  values  than  it  did  for  the  elliptic  leading  edge.  In  addition,  the 
pressure  plateau  is  much  longer  and  higher  for  the  circular  leading  edge.  This  is,  again, 
because  of  the  greater  deceleration  that  occurs  around  the  circular  leading  edge,  which 
resulted  in  the  fonnation  of  a  separation  bubble.  With  an  elliptic  leading  edge,  a 
separation  bubble  may  not  develop  at  zero  incidence  because  separation  is  avoided 
entirely.  When  the  bubbles  do  fonn,  they  are  usually  shorter  in  length  and  are  able  to 
reattach  to  the  blade  surface  at  a  higher  incidence  than  those  which  formed  on  the  circular 
leading-edge  models.  The  boundary  layer  momentum  thickness  on  the  circular  leading- 
edge  model  was  significantly  larger  than  that  of  the  elliptic  leading  edge,  and  indicates 
that  an  elliptic  leading  edge  reduces  losses. 

Many  airfoil  sections  used  in  compressors  and  turbines  have  discontinuity  in 
curvature  of  the  leading  edge  [1]  at  the  transition  point  from  leading  edge  to  airfoil  body. 
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At  the  very  front  of  the  blade,  the  curvature  is  quite  small,  and  drastically  increases  at  the 
point  where  the  leading  edge  joins  the  body  of  the  blade.  Typical  compressor  blades,  such 
as  the  NACA  65  series,  have  a  maximum  thickness  that  ranges  from  roughly  three  to  12 
percent  of  the  chord  length,  and  a  leading-edge  thickness  that  is  around  0.4  percent  of  the 
chord  [1].  When  compared  to  the  radius  of  curvature  of  the  leading  edge,  the  radius  of  the 
blade  body  is  almost  flat.  This  discontinuity  in  curvature,  such  as  with  the  circular 
leading-edge  model,  results  in  a  separation  bubble  forming  even  at  zero  incidence. 

Huoxing  et  al.  [2]  also  determined  that  the  leading-edge  shape  has  a  large 
influence  on  flow  characteristics.  This  affects  separation  and  transition  of  the  flow,  as 
well  as  the  boundary  layer  properties  after  the  reattachment  point.  In  their  experiments, 
they  tested  a  circular  leading-edge  profde,  as  well  as  a  model  with  a  leading-edge 
ellipticity  ratio  of  2.5.  In  tests  of  their  elliptic  leading-edge  model,  the  suction  peak  in  the 
pressure  coefficient  distribution  was  much  lower,  and  the  flow  did  not  separate  even  at 
six  degrees  angle  of  attack.  They  also  witnessed  that  the  effect  of  the  leading-edge  shape 
on  the  pressure  side  of  the  blade  was  rather  weak  when  compared  to  the  influence  on  the 
flow  over  the  suction  side.  Even  at  low  Reynolds  numbers,  Huoxing  et  al.  [2]  found  that  a 
separation  bubble  could  form  at  low  incidence  on  the  circular  leading-edge  model. 
Comparatively,  if  a  bubble  did  form  on  the  elliptic  leading-edge  model,  it  was  shorter, 
and  allowed  the  flow  to  reattach  at  higher  incidences  than  was  possible  with  the  circular 
leading  edge.  Huoxing  et  al.  [2]  define  a  separation  bubble  as  a  region  where  a  strong 
adverse  pressure  gradient  can  cause  a  laminar  boundary  layer  to  separate.  However,  once 
the  flow  is  separated,  the  shear  layer  becomes  highly  unstable  and  can  undergo  rapid 
transition  to  turbulent  flow.  This  turbulent  flow  in  the  separation  bubble  is  able  to  entrain 
more  fluid  from  surrounding  laminar  flow  and  the  turbulent  shear  layer  then  usually 
bends  back  towards  the  surface  of  the  blade  and  reattaches  as  a  turbulent  boundary  layer. 
Their  tests  also  revealed  that  the  length  of  the  separation  bubble,  as  well  as  the  pressure 
distribution  and  boundary  layer  characteristics  downstream  of  the  bubble  were  strong 
functions  of  incidence.  As  the  flow  reattached  to  the  surface  there  was  a  steep  pressure 
rise.  As  they  increased  testing  angle  of  attack,  the  formation  of  the  bubble  occurred  at  the 
same  chord-wise  point,  however  the  bubble  became  longer  and  the  level  of  uniformity  in 
the  pressure  plateau  became  lower. 
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Chandrasekhara  et  al.  [3]  perfonned  experiments  on  a  VR-12  helicopter  rotor 
blade  during  which  they  statically  and  dynamically  drooped  the  first  25  percent  of  the 
airfoil  from  zero  degrees  up  through  20  degrees.  Their  experiments  were  performed  at 
relatively  low  Mach  numbers,  from  0.1  to  0.4.  The  droop  angle  is  denoted  as  8 .  They 
discovered  that  drooping  the  leading  edge  allowed  for  better  control  of  the  flow  around 
the  nose.  When  they  compared  the  pressure  distribution  over  the  drooped  airfoil  to  that  of 
the  undrooped  airfoil,  the  pressure  peak  was  dramatically  different.  On  the  undrooped 
blade,  the  pressure  coefficient  peaked  out  at  -7.92,  where  as  on  one  of  their  statically 
drooped  models,  it  only  reached  -3.46  for  the  same  testing  conditions.  Also,  their  study 
determined  that  the  suction  side  pressure  values  of  a  drooped  leading-edge  model 
appeared  to  lag  the  basic  airfoil  in  angle  of  attack.  Additionally,  when  the  leading  edge 
was  drooped,  the  flow  remained  attached  to  a  point  farther  down  the  blade.  The  research 
conducted  by  Chandrasekhara  et  al.  [3]  also  illustrates  the  importance  of  a  smooth  blend 
between  the  drooped  leading  edge  and  the  rest  of  the  airfoil  section.  For  example,  they 
observed  a  slight  distortion  in  the  pressure  distribution  near  the  blend  point  that  resulted 
from  the  small  backward  facing  step  present  where  the  hinge  connected  the  drooped 
leading  edge  to  the  main  body.  This  motivates  special  attention  to  be  paid  to  blending  the 
leading  edge  to  avoid  tripping  the  flow. 

The  benefits  of  a  drooped  leading  edge  are  particularly  noticeable  at  higher  angles 
of  attack.  Chandrasekhara  et  al.  [3]  were  able  to  achieve  fully  attached  leading-edge  flow 
even  at  20  degrees  angle  of  attack  with  a  drooped  nose.  Drooping  the  leading  edge  also 
allows  for  greater  control  of  the  acceleration  of  the  flow,  similar  to  the  effects  of  an 
elliptic  edge  when  compared  to  that  of  a  circular  leading  edge.  Their  drooped  model  was 
able  to  keep  the  flow  well  subsonic  when  compared  to  the  undrooped  model,  even  at  20 
degrees  alpha.  The  peak  suction  over  the  basic  airfoil  rose  steeply  with  angle  of  attack 
[3].  The  flatter  pressure  distribution  seen  with  a  drooped  leading-edge  implied  that  there 
is  still  a  reasonably  high  amount  of  lift  produced.  Drooping  the  leading  edge  was  able  to 
largely  mitigate  compressibility  affects,  which  would  enable  a  better  performance  to  be 
attained  from  the  same  airfoil  at  higher  Mach  numbers.  Drooping  the  leading  edge  also 
pushes  the  angle  of  attack  at  which  the  VR-12  airfoil  stalls  to  a  higher  value  than  attained 
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with  the  undrooped  leading  edge.  The  data  indicated  that  a  greater  operating  envelope  is 
attainable  with  a  drooped  leading  edge. 

Tuck  and  Dostovalova  [4]  found  that  drooping  the  leading  edge  by  sharpening  the 
nose  can  result  in  better  flow  control.  In  their  studies,  they  changed  the  shape  of  the  nose 
to  best  delay  separation  of  the  laminar  boundary  layer.  Park  [5]  arrived  at  the  same 
conclusion  and  asserts  that  a  wedge  shaped  leading  edge  can  be  inexpensively  created 
using  casting  technology  that  leaves  a  relatively  thick  leading  edge.  Afterwards,  the  blade 
can  easily  be  retouched  to  create  the  desired  leading-edge  shape  using  milling. 

Garzon  and  Darmafol  [6]  conducted  studies  on  variations  in  leading-edge 
geometry  that  might  arise  due  to  manufacturing  inaccuracies.  The  greatest  losses 
occurred  as  a  result  of  variations  in  the  first  five  percent  of  the  blade,  illustrating  the 
dramatic  effects  that  the  leading-edge  shape  and  direction  can  have  on  determining  flow 
characteristics  over  the  rest  of  the  airfoil.  Additionally,  they  observed  substantial 
increases  in  boundary  layer  momentum  thickness  as  a  result  of  these  geometric 
variations. 

Tuck  [7]  also  performed  a  study  that  compared  lift  results  to  those  listed  in  Abbott 
and  Von  Doenhoff  [8],  and  found  that  the  reattachment  of  the  flow  is  responsible  for  the 
underprediction  of  the  angle  of  attack  for  maximum  lift  by  a  factor  of  0.5  to  0.8.  He  also 
asserts  that  the  shape  of  the  leading  edge  will  greatly  affect  whether  or  not  the  flow 
reattaches.  Tuck  also  found  that  sharp-edged  wings  stall  very  quickly,  further  illustrating 
the  impact  of  leading-edge  geometry  on  attaining  smooth,  attached  flow  over  the  airfoil. 

Suder  et  al.  [9]  found  that  increasing  the  thickness  and  adding  roughness  to  a 
transonic  axial-flow  compressor  rotor  can  result  in  a  significant  degradation  in  overall 
performance.  In  manufacturing  and  repairing  transonic  axial  blades,  Suder  et  al.  [9] 
believe  that  the  goal  should  be  to  make  the  leading  edge  as  smooth  as  possible  and  so  that 
it  will  stay  that  way  for  as  long  an  operating  period  as  possible.  Applying  a  protective 
coating  to  the  blade  increases  resistance  to  degradation.  These  coatings  do,  however, 
results  in  an  increase  of  the  suction  side  boundary  layer  thickness  due  to  their  impact  on 
leading-edge  geometry. 
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The  unintentional  blunting  of  the  leading  edge  can  have  detrimental  effects  on  the 
flow  field  over  an  airfoil.  Reid  and  Urasek  [10]  found  that  for  relatively  small 
compressors,  the  leading  edge  of  the  rotor  blades  may  be  blunt  due  to  manufacturing 
tolerances.  For  larger  blades,  the  leading  edge  may  become  blunt  due  to  foreign  object 
damage  or  corrosion.  In  transonic  rotors,  a  blunt  leading  edge  tends  to  increase  the  shock 
losses  associated  with  a  detached  bow  wave  system.  With  a  blunt  leading  edge,  the  total 
pressure  ratio,  total  temperature  ratio  and  adiabatic  efficiency  is  much  lower  than  that  for 
the  circular  leading  edge. 

Tain  and  Cumpsty  [11]  performed  experiments  and  simulations  with  a  circular 
leading  edge  attached  to  a  flat  airfoil  at  Mach  numbers  ranging  from  0.6  to  0.95.  The 
extent  to  which  the  airflow  accelerates  around  the  leading  edge  is  dependent  on  the  shape 
of  the  leading  edge.  They  found  that  you  can  expect  the  flow  to  become  sonic  around  a 
circular  leading  edge  of  an  airfoil. 

Computational  Fluid  Dynamics  is  heavily  relied  upon  as  a  method  of  predicting 
leading-edge  flow  behavior.  Calvert  [12]  investigated  an  inviscid-viscous  method  to 
modeling  leading-edge  separation  bubbles,  and  found  that  they  are  most  likely  to  occur 
near  the  leading  edge  of  sharp-edged  blade  sections  in  axial  compressors  and  turbines. 
This  is  particularly  the  case  when  the  blade  sections  are  operated  at  positive  incidence. 
Calvert  [12]  also  found  that  the  overall  performance  of  the  blade  section  can  be 
significantly  affected  by  the  thickness  of  the  initial  boundary  layer. 

Mason  [13]  asserts  that  the  Navier-Stokes  solutions  used  in  computational  fluid 
dynamic  studies  do  not  adequately  address  the  flow  field  solution  near  the  leading  edge. 
His  work  stresses  the  importance  of  leading-edge  flow  on  boundary  layer  development. 

Environmental  conditions  have  been  shown  to  significantly  alter  the  leading-edge 
shape.  Wickens  and  Nguyen  [14]  observed  leading-edge  changes  that  resulted  from  ice 
accretion  on  a  propeller,  and  found  that  there  was  a  large  loss  of  lift  and  an  increase  in 
form  drag  as  a  result  of  the  modification  of  the  leading  edge.  Their  tests  also  showed  a 
large  dependency  on  Reynolds  number. 

Simulations  of  ice  accretion  were  performed  by  Kwon  and  Sankar  [15]  using  a 
three-dimensional  compressible  Navier-Stokes  solver.  Their  investigation  indicates  that 
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leading-edge  separation  and  reattachment  can  be  caused  by  changes  to  the  leading-edge 
shape  from  ice  accretion.  Similarly,  Broeren  et  al.  [16]  confirmed  that  the  separation 
bubble  became  larger  as  the  leading-edge  shape  became  increasingly  non-uniform  as  a 
result  of  ice  accretion. 

Lu  and  Xu  [17]  performed  an  investigation  of  a  circular  leading  edge  with  a  flat 
section  both  experimentally  and  using  CFD.  They  found  that  the  double  suction  spikes  on 
the  flat  edges  are  much  weaker  than  the  single  spike  that  appears  with  a  traditional 
circular  leading  edge.  At  higher  incidences,  the  separation  downstream  of  the  leading 
edge  was  weaker  with  the  flat.  Lu  and  Xu  [17]  claim  that  this  improvement  in 
performance  was  better  or  equal  to  that  expected  with  an  elliptical  leading  edge;  however 
they  believe  that  the  circular  leading  edge  with  a  flat  is  more  tolerant  to  manufacturing 
inaccuracy  than  the  elliptical  shaped  leading  edge. 

2.  Summary  of  Literature  Survey 

Several  researchers  have  experimentally  and  numerically  investigated  the  effects 
of  leading-edge  geometry  modification  on  flow  field  properties.  They  have  concluded 
that  ellipticity  and  droop  angle  can  significantly  impact  flow  separation  and  transition  to 
turbulence.  The  combined  effects  on  pressure  distribution  and  aerodynamic  losses  of  a 
drooped  elliptic  leading  edge  at  varying  incidences  have  not  been  extensively  studied  in 
the  flow  conditions  present  in  a  typical  compressor. 


C.  OBJECTIVES 

The  current  study  focuses  on  the  combined  effects  of  ellipticity  ratio  and  leading- 
edge  droop  over  a  range  of  incidences  on  the  aerodynamic  losses  for  a  compressor  blade. 
It  is  the  overall  aim  to  provide  insight  into  how  these  geometrical  variations  impact  the 
flow  at  a  Mach  number  and  Reynolds  number  typical  of  those  in  a  compressor.  More 
specifically,  the  objectives  of  this  thesis  are  to: 

1)  Demonstrate  ESI-CFD’s  ability  to  accurately  predict  flow  behavior  when 
compared  to  results  of  previous  experimental  studies. 
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2)  Obtain  numerical  results  for  the  pressure  distribution  over  the  blade  as  well  as 
boundary  layer  momentum  thickness  at  a  specified  point  on  the  model  for  simulations 
performed  at  typical  compressor  conditions. 

3)  Observe  the  impact  of  incidence  on  the  flow  over  models  with  leading-edge 
ellipticity  and  droop. 

4)  Determine  an  optimum  leading-edge  ellipticity  ratio  and  droop  angle  that 
provides  the  greatest  reduction  in  aerodynamic  losses. 


D.  ORGANIZATION 

Chapter  II  presents  the  methodology  used  to  conduct  the  research  on  leading-edge 
shape  and  direction.  Methods  used  in  model  creation  are  discussed.  Also,  the  steps 
involved  in  the  CFD  simulations,  as  well  as  the  simulation  conditions  and  variations  are 
explained. 

Chapter  III  provides  a  graphical  representation  of  the  results  obtained  during  the 
simulations.  A  discussion  is  also  presented  to  further  explain  and  identify  key  results  of 
the  study  and  their  significance. 

Chapter  IV  contains  conclusions  that  are  drawn  based  on  the  results  obtained,  and 
a  summary  of  their  impact  on  the  study’s  objectives  is  also  provided. 
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II.  METHODOLOGY 


A.  COMPUTATIONAL  FLUID  DYNAMICS 

The  research  conducted  in  the  current  study  was  performed  using  Computation 
Fluid  Dynamics.  The  ESI-CFD  program,  developed  by  CFDRC,  was  chosen  for  use  in 
the  simulations.  This  code  is  based  on  three  main  parts;  a  pre-processor,  the  processor 
and  the  post-processor.  These  subprograms  were  used  to  define  the  model,  setup  and  run 
the  problem,  and  collect  the  data  that  the  simulations  produced. 

The  model  geometries  were  created  in  MATLAB  and  coordinates  for  plates  and 
their  leading  edges  were  inputted  into  the  pre -processor,  CFD-GEOM.  In  this 
subprogram,  geometrical  parameters  of  the  simulation  were  established  such  as  where  the 
inlet  would  be  in  relation  to  the  model.  After  creating  the  boundaries  of  the  problem,  a 
grid  was  be  generated  which  allows  the  user  to  discretize  the  geometry  into 
computational  cells. 

In  the  processor,  CFD-ACE+,  initial,  volume,  and  boundary  conditions  were 
applied  to  the  model  created  in  CFD-GEOM.  Flow  conditions  were  defined  and 
preferences  were  set  for  spatial  differencing,  iteration  limits,  as  well  as  relaxation  factors. 
The  problem  was  run  from  within  CFD-ACE+,  and  allowed  the  user  to  observe  how  well 
the  problem  converges. 

The  post  processor,  CFD-VIEW,  allows  the  user  to  visualize  and  analyze  the 
solution.  Surface  and  line  plots  were  created  to  retrieve  data  from  within  the  flow  field. 
The  data  sought  by  the  user  was  then  taken  from  the  post-processor  and  reduced  in 
MATLAB. 

B.  VALIDATION  OF  CFD-ACE+ 

To  validate  the  results  obtained  using  the  ESI-CFD  programs,  a  comparison  study 
was  performed  based  on  the  work  of  Walraevens  and  Cumpsty  [1],  A  simulation  model 
was  constructed  and  run  at  the  same  conditions  as  an  experiment  that  they  perfonned  on 
an  elliptic  leading-edge  plate,  as  outlined  in  Reference  [1].  Table  1  lists  the  conditions  at 
which  the  simulation  and  experiment  were  run. 
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Table  1.  Conditions  For  Walraevens-Cumpsty  Test 


Ellipticity  Ratio 

1.89 

Angle  of  Attack 

0  deg 

Velocity 

7.2227  m/s 

Density 

1.1614  kg/mA3 

Kinematic  Viscosity 

1.589E-5  mA2/s 

Reference  Static  Pressure 

0.1  MPa 

Pressure  data  were  taken  around  the  surface  of  the  model.  This  was  used  in 
obtaining  the  pressure  coefficient.  Figure  8  compares  the  pressure  coefficients  obtained 
by  Walraevens  and  Cumpsty  to  that  obtained  by  CFD. 
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Surface  Distance  From  Leading  Edge,  meters 

Figure  8.  Comparison  of  Walraevens-Cumpsty  and  CFD  Pressure  Coefficients  (-Cp)  for  a 

Flat  Plate  with  Leading-Edge  Ellipticity  Ratio  of  1.89 

As  Figure  8  shows,  the  pressure  plots  do  not  exactly  coincide.  The  disparity 

between  the  two  pressure  plots  lies  in  the  fact  that  in  the  experimental  setup,  Walraevens 

and  Cumpsty  used  0.5  percent  turbulence;  however,  in  the  CFD  simulation,  no  turbulence 
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was  used.  However,  these  results  are  able  to  validate  the  CFD  code’s  prediction  ability 
because  the  pressure  peak  has  approximately  the  same  value  between  the  simulation  and 
the  experiment.  Also,  the  location  of  the  pressure  peak  along  the  surface  of  the  plate  is 
almost  identical  in  the  experiment  and  simulation.  This  suggests  that  the  CFD  simulations 
provide  results  that  are  close  to,  if  not  the  same  as,  what  would  be  obtained 
experimentally. 

C.  MODEL  CREATION 

1.  Ellipticity  and  Droop 


c 


ER=I 


ER=5.5 


Figure  9.  Flat  Plate  Models  with  Leading  Edges  of  Varying  Ellipticity  Ratio 
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Figure  9  shows  the  physical  impact  of  varying  the  ellipticity  ratio  of  the  leading 
edge.  The  models  shown  in  the  figure  have,  from  top  to  bottom,  an  ellipticity  ratio  of  1, 
2.5,  4,  and  5.5.  The  calculation  of  the  circumference  of  the  ellipse  for  the  undrooped 
models  was  done  using  Ramanujan’s  approximation  [18].  The  equation  used  was: 

C  =  x([3a  +  3b-j(a  +  3b)(b  +  3a)]j  (1) 


Equation  (1)  allowed  for  the  determination  of  the  distance  along  the  surface  from 
the  geometric  leading  edge  to  the  blend  point  of  the  undrooped  models.  Figure  49  of 
Appendix  E  presents  this  equation  graphically  for  the  undrooped  models  tested. 


0  degrees  droop 


3.33  degrees  droop 


6.66  degrees  droop 


10.0  degrees  droop 


13.3  degrees  droop 


Figure  10.  Flat  Plate  Models  with  Varying  Droop  Angles  of  0,  3.33,  6.66,  10.0,  and  13.3 

degrees 
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The  effect  of  droop  is  shown  in  Figure  10.  Each  blade  shown  has  the  same 
ellipticity  ratio;  however  the  droop  varies  from,  top  to  bottom,  0,  3.3,  6.6,  10,  and  13.3 
degrees.  The  droop  angle  is  defined  as  the  angle  from  the  centroid  of  the  ellipse  to  the 
apex  of  the  semimajor  axis  at  the  leading  edge. 

Each  model  was  created  by  taking  a  flat  plate  of  constant  thickness  and  attaching 
a  different  leading  edge  to  it.  Although  not  a  full  representation  of  the  geometry  of  an 
entire  compressor  blade,  using  a  flat  plate  as  the  body  still  allows  for  the  analysis  of  the 
leading  edge’s  effects  on  the  flow  characteristics  over  the  model,  and  what  could  be 
expected  for  a  compressor  blade.  The  circular  leading-edge  model  required  that  a  circle 
of  the  same  diameter  as  the  plate  be  used,  and  simply  attached  as  a  half  circle  to  the  front 
of  the  blade  model.  A  similar  method  was  used  for  the  undrooped  elliptic  leading-edge 
models.  However,  for  the  models  where  the  leading  edge  was  drooped,  further 
modifications  were  required.  The  thickness  of  the  plate  body  is  a  constant  eight 
millimeters,  and  the  length  of  the  body,  without  the  leading  edge,  is  0.0945  meters. 

2.  Blending  of  the  Leading  Edge 

The  drooped  leading  edges  were  created  by  taking  an  ellipse  of  specified 
ellipticity  ratio,  and  rotating  it  to  the  desired  droop  angle.  This  canted  ellipse  was  then 
laid  over  the  front  of  the  flat  plate.  Initially,  the  excess  ellipse  and  flat  plate  that  remained 
were  “trimmed”  away,  and  a  complete  airfoil  was  created. 
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Figure  1 1 .  Drooped  Leading  Edge  Before  Blending 

As  Figure  1 1  shows,  the  intersection  of  the  drooped,  elliptic  leading  edge  and  the 
flat  plate  can  clearly  be  seen.  The  resulting  combination  of  the  two  geometries  creates  a 
sharp  point  on  both  the  suction  and  pressure  sides.  These  pronounced  comers  can  trip  the 
boundary  layer  flow,  and  cause  early  transition  to  turbulence  or  the  formation  of  a 
separation  bubble. 

In  the  design  of  leading-edge  shapes,  a  very  important  element  to  take  into 
account  is  the  radius  of  curvature  of  the  leading  edge,  as  well  as  how  it  changes  from  the 
leading  edge  towards  the  body.  A  circular  leading  edge  has  a  well  defined,  constant 
radius.  At  the  point  where  the  circular  leading  edge  is  mated  with  the  airfoil  body,  there  is 
a  very  large  increase  in  radius  of  curvature.  This  dramatic  change  in  curvature  has  been 
found  to  cause  separation.  However,  with  an  elliptical  leading  edge,  the  rate  at  which  the 
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radius  of  curvature  changes  is  not  as  dramatic.  Figure  12  presents  a  comparison  of  the 
radius  of  curvature  at  the  leading  edge  and  the  blend  point  for  leading  edges  of  varying 
elliptic ity  ratio. 


Figure  12.  Variation  of  Radius  of  Curvature  with  Arclength  From  Geometric  Leading  Edge 

to  Blend  Point  at  Semiminor  Apex 

At  the  very  front  of  an  elliptic  leading  edge  the  radius  is  small,  and  gradually 
increases  towards  the  apex  of  the  semiminor  axis.  In  Figure  11,  where  there  is  a  sharp 
comer,  a  very  small  radius  of  curvature  exists.  In  order  to  allow  for  a  more  gradual 
change  in  radius  of  curvature,  the  sharp  corner  where  the  leading  edge  and  blade  body  are 
joined  was  replaced.  This  was  achieved  by  removing  a  section  of  the  leading  edge  and 
blade  body  where  the  two  join  on  both  the  suction  and  pressure  sides.  Around  this  point, 
a  section  of  a  circular  arc  was  fitted.  As  a  result,  a  smooth  transition  is  created  from  the 
leading  edge  to  the  body  of  the  blade,  and  the  likelihood  that  the  transition  point  between 
the  two  would  result  in  early  transition  to  turbulence  or  a  separation  bubble  is  decreased. 
Figure  13  illustrates  the  newly  fitted  transition  region  from  the  leading  edge  to  the  blade 
body.  The  two  lines  running  vertically  join  at  the  center  of  the  circular  section  that  was 
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used  to  create  the  transition  on  the  pressure  side  of  the  blade,  and  indicates  the  typical 
radius  of  curvature  used  for  these  transition  sections. 


Figure  13.  Radius  of  Curvature  for  Pressure  Side  Blending  Region  of  Drooped  Elliptic 

Leading-Edge  Attached  to  a  Flat  Plate 

Figure  14  shows  a  close-up  comparison  of  the  sharp  leading  edge  originally  used 
to  the  circular  arc  fitted  to  the  models  for  a  smooth  transition.  Figure  15  shows  an 
example  of  a  final  leading-edge  shape  that  was  used  after  being  refitted  with  a  smoother 
joining  region.  Instead  of  the  sharp  corners  that  are  shown  in  Figure  11,  the  leading  edge 
is  smoothly  mated  to  the  airfoil  body.  As  the  droop  angle  increased,  smaller  radii  of 
curvature  sections  were  fitted.  Also,  for  the  larger  ellipticity  ratio  blades,  larger  radii  of 
curvature  sections  were  used  because  of  the  larger  radius  of  curvature  that  occurs  at  the 
semiminor  axis  with  a  larger  ellipse. 
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Figure  14.  Comparison  of  Original  Leading-Edge  Transition  to  Blended  Leading-Edge 

Transition  on  Pressure  Side  of  Plate 


Figure  15.  Final  Blade  After  Blending  of  Leading  Edge  on  Pressure  and  Suction  Sides 

The  schematic  in  Figure  15  depicts  the  final  product  for  a  leading  edge  with  an 
ellipticity  ratio  of  2.5  and  a  droop  angle  of  13.3  degrees.  The  transition  regions  on  both 
the  suction  and  pressure  side  are  blended  to  allow  for  the  smooth  flow  from  the  leading 
edge  to  the  plate  body. 
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Table  2  lists  the  leading-edge  geometric  variations  that  were  studied.  Ellipticity 
ratios  were  first  chosen,  ranging  in  value  from  less  than  one  to  5.5.  Perfonning 
simulations  on  a  blade  model  with  an  ellipticity  ratio  of  less  than  one  allows  for  the 
effects  of  a  blunt  leading  edge  to  be  observed.  Past  theory  holds  that  an  ellipticity  ratio  of 
four  is  optimum  for  the  minimization  of  losses.  Therefore,  models  with  an  ellipticity  ratio 
lower  and  higher  than  this  were  tested  to  determine  if  an  ellipticity  ratio  of  four  is  indeed 
optimal. 


Tal 


fie  2.  Model  Ellipticity  Ratios  and  Droop  Angles 


Ellipticity  Ratio 

Droop  Angle  (degrees) 

1 

0.0 

2.5 

0.0 

2.5 

3.33 

2.5 

6.66 

2.5 

10.0 

2.5 

13.3 

3 

0.0 

3.5 

0.0 

4 

0.0 

4 

3.33 

4 

6.66 

4 

10.0 

4 

13.3 

5.5 

0.0 

5.5 

3.33 

5.5 

6.66 

5.5 

10.0 

5.5 

13.3 

After  the  geometry  of  the  models  to  be  used  had  been  set,  simulation  models  had 
to  be  constructed  using  the  pre-processor,  CFD-GEOM.  ‘Tab’  files  were  made  to  input 
the  model  geometry  data  into  this  subprogram.  Edges  and  surfaces  were  defined,  and  a 
grid  was  constructed  based  on  the  basic  plate  and  leading-edge  geometry. 
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Figure  16.  Sample  Grid  Generated  in  CFD-GEOM  Used  For  Simulations 


Figure  16  represents  a  typical  grid  used  to  define  the  simulation  space  about  the 
plate  and  leading  edge.  Despite  a  variation  in  overall  blade  length  due  to  different 
ellipticity  ratios,  the  grid  size  remained  constant.  The  overall  height  and  length  of  the  grid 
are  0.2  meters  and  0.2995  respectively.  The  inlet  region  is  semicircular,  and  was  designed 
this  way  so  that  the  grid  contours  would  more  closely  match  the  general  leading-edge 
shape.  On  both  the  suction  side  and  pressure  side  of  the  plate,  140  grid  points  were  used, 
and  80  grid  points  placed  were  placed  along  the  leading  edge.  Blade  length  was  held  at  a 
constant  0.0945  meters  between  the  different  models,  and  the  length  of  the  leading  edge 
was  varied  according  to  the  specified  ellipticity  ratio.  Defining  the  boundaries  with  the 
number  of  grid  points  specified  in  Figure  16  resulted  in  35,800  nodes  used  for 
simulations. 

Figure  17  provides  a  close  up  view  of  the  grid  at  the  leading  edge  of  a  model.  The 
grid  becomes  much  denser  close  to  the  surface  of  the  blade,  increasing  with  a  power  law 
factor  of  four.  Grid  refinement  was  necessary  close  to  the  blade  surface  to  allow  to  high 
spatial  resolution  within  the  boundary  layer.  To  ensure  that  a  sufficient  number  of  grid 
points  would  exist  within  the  boundary  layer,  the  boundary  layer  momentum  thickness 
was  calculated  for  a  flat  plate  of  the  same  length  as  the  blade  body.  At  the  point  where 
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boundary  layer  momentum  thickness  was  to  measured,  a  minimum  of  ten  grid  points 
were  contained  within  the  boundary  layer  for  all  cases. 


Figure  17.  Close-up  of  Grid  At  Leading  Edge 


A  grid  dependency  study  was  conducted  to  determine  the  impact  of  the  grid 
resolution  used  for  the  simulations.  This  was  perfonned  by  doubling  the  number  of  grid 
points  on  each  boundary  and  comparing  results  for  boundary  layer  momentum  thickness. 
Although  the  answers  between  the  two  grids  were  different,  their  disparity  was  less  than  a 
few  percent,  validating  the  original  grid’s  ability  to  produce  reliable  results  without  the 
dramatic  increase  in  computation  time. 

D.  SIMULATION  CONDITIONS 

Before  simulation  runs  could  be  made  using  CFD-ACE+,  the  problem  had  to  be 
properly  set  up  with  defined  boundary  conditions,  volume  conditions,  and  solver 
conditions.  Table  3  lists  the  conditions  that  were  inputted  into  CFD-ACE+.  A  plate  length 
Reynolds  number  of  about  1 .4  million  was  detennined  according  to  the  equation 
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(2) 


A  Mach  number  of  0.65  is  typical  of  core  compressor  conditions.  The  Mach 
number  was  arrived  at  according  to  equation  (3) 


Ma  = 


(3) 


The  kinematic  viscosity  listed  was  held  constant  for  all  simulations.  Boundary 
conditions  were  defined  for  the  inlet  and  outlet.  The  top  and  bottom  boundaries  in  Figure 
12  were  set  to  be  periodic.  This  allowed  the  solution  for  the  simulation  of  one  blade  to 
represent  what  could  be  expected  for  a  cascade  of  compressor  blades.  The  inlet  condition 
was  set  according  to  the  velocities  listed  in  Table  4,  to  impact  the  specified  incidence  on 
the  model.  The  initial  condition  was  set  to  be  close  to  the  inlet  velocity  to  allow  for  faster 
convergence.  Central  spatial  differencing  was  used  with  a  blending  factor  of  0.1. 
Relaxations  were  set  to  2  for  velocities  and  0.1  for  pressure.  A  convergence  criterion  was 
set  to  0.00001.  All  blade  models  were  run  at  incidences  of  0,  3,  6,  10  and  13  degrees. 


Table  3.  Simulation  Test  Conditions 


Temperature 

293  K 

Density 

1.2047  kg/mA3 

Kinematic  Viscosity 

1.508E-5  mA2/s 

Reference  Static  Pressure 

101325  Pa 

Table  4.  Simulation  Component  Velocities 


Incidence  (deg) 

U-Component  Velocity  (m/s) 

V-Component  Velocity  (m/s) 

0 

223.02 

0 

3 

222.72 

11.67 

6 

221.80 

23.31 

10 

219.64 

38.73 

13 

217.31 

50.17 
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III.  RESULTS  AND  DISCUSSION 


A.  RESULTS 

1.  Pressure  Distribution 

The  pressure  coefficient  was  determined  for  each  run  by  taking  pressure  data 
around  the  surface  of  the  blade,  and  reduced  using  the  equation  (4) 


C, 


P-P„ 


(4) 


This  equation  is  most  commonly  used  for  low  speed,  incompressible  flows.  The 
denominator  of  equation  (4)  is  often  represented  by  the  difference  between  the  total  and 
static  pressure,  Pn-P,  which  is  approximately  equal  to  the  dynamic  head  for  low  Mach 
number  flows.  However,  because  this  study  used  a  consistent  free  stream  Mach  number 
of  0.65,  this  equation  provides  very  close  answers  even  in  the  compressible  flow  region. 

A  compilation  of  pressure  distributions  is  included  in  this  section.  It  is  a  selection 
of  figures  that  represent  the  general  behavior  of  the  pressure  coefficient  over  the  models. 
Plots  for  each  individual  model  can  be  found  in  Appendix  C.  The  figures  shown  represent 
the  pressure  coefficient  near  the  leading  edge  of  the  blade,  where  the  flow  characteristics 
have  a  significant  effect  on  the  flow  over  the  rest  of  the  blade.  Measurements  are  based 
on  the  distance  from  the  nose  of  the  leading  edge. 
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Figure  18.  Pressure  Distribution  (-Cp)  Over  Circular  Leading-Edge  Model  at  Various  Angles 

of  Attack  with  Blendpoint  Shown 

One  of  the  goals  of  this  study  was  to  determine  a  leading-edge  geometry  that 
provided  a  better  pressure  distribution  over  the  blade.  It  was  therefore  necessary  to 
analyze  the  pressure  distribution  over  a  model  with  a  circular  leading  edge,  such  as  is 
presented  in  Figure  18.  The  figure  shows  that  as  the  angle  of  attack  increases,  the  height 
of  the  pressure  spike  increases,  and  the  flow  farther  downstream  of  the  leading  edge 
becomes  increasingly  poorly  behaved.  If  the  plate  were  infinitely  long,  the  pressure 
coefficient  should  eventually  go  to  zero.  Also,  from  the  line  representing  the  pressure 
distribution  at  zero  degrees  angle  of  attack,  obvious  perturbations  in  the  pressure 
coefficient  can  be  seen.  This  indicates  a  strong  adverse  pressure  gradient  exists  on  the 
circular  leading  edge,  and  that  the  flow  may  be  separating  even  at  zero  incidence. 
Visualizations  of  the  flow  from  the  post-processor  confirmed  this,  as  well  as  the  results 
obtained  by  Walraevens  and  Cumpsty  [1]  in  their  experiments  using  a  circular  leading 
edge  attached  to  a  flat  plate. 
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Surface  Distance  from  Geometric  Leading  Edge,  meters 


Figure  19.  Pressure  Distribution  (-Cp)  Over  All  Ellipticities  with  Zero  Droop  at  Zero 

Degrees  Incidence  with  Blendpoints  Shown 

Figure  19  displays  the  pressure  variations  as  a  function  of  ellipticity  ratio.  The 
important  trend  to  notice  in  this  figure  is  that  the  pressure  spike  decreases  its  magnitude 
as  the  ellipticity  ratio  increases.  Also,  the  region  of  pressure  elevation  becomes  much 
more  elongated  with  the  more  elliptic  leading  edges.  This  is  caused  by  the  more  gradual 
acceleration  and  deceleration  of  the  flow  around  the  leading  edge  as  the  ellipticity 
increases.  The  backside  of  the  pressure  peak  is  also  important,  and  is  a  strong  indication 
of  whether  or  not  the  flow  separates.  The  plateaus  that  appear  on  the  backside  of  the 
pressure  peak  on  the  circular  leading  edge  suggest  the  presence  of  a  separation  bubble. 
Even  for  the  model  with  an  ellipticity  ratio  of  2.5,  a  small  perturbation  and  plateau  can  be 
seen  on  the  backside  of  the  pressure  peak.  This  indication  of  flow  separation  was 
confirmed  with  flow  visualization  in  the  post-processor.  However,  for  the  blade  model 
with  an  ellipticity  ratio  of  5.5,  for  example,  the  backside  of  the  pressure  peak  is  very 
smooth,  and  suggests  that  the  flow  remains  fully  attached  to  the  blade  surface  and  does 
not  separate. 
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Figure  20.  Pressure  Distribution  (-Cp)  Over  Model  with  an  Ellipticity  Ratio  of  2.5  with  No 

Droop  at  Various  Angles  of  Attack  with  Blendpoint  Shown 

Observing  the  impact  of  ellipticity  on  pressure  distribution,  Figure  20  provides 
insight  into  the  flow  behavior  for  the  blade  model  with  an  ellipticity  ratio  of  2.5.  Similar 
to  the  trend  observed  in  Figure  18  for  the  circular  leading  edge  model,  as  the  angle  of 
attack  increases,  the  height  of  the  pressure  spike  increases.  Additionally,  the  flow  on  the 
suction  side  is  well  behaved  at  the  lower  incidences,  as  indicated  by  the  smooth  pressure 
distribution  after  the  peak  pressure  point.  Despite  the  fact  that  at  the  higher  incidences 
Figure  20  suggests  separation,  the  pressure  peak  is  still  significantly  lower  than  it  was  for 
the  circular  leading-edge  model.  The  oscillations  after  the  pressure  peak,  as  shown  for  the 
higher  incidences,  indicate  that  there  are  multiple  regions  of  separation  on  the  suction 
side  of  the  plate. 
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Figure  21.  Pressure  Distribution  (-Cp)  Over  Model  with  an  Ellipticity  Ratio  of  5.5  and  Ten 

Degrees  Droop  at  Various  Angles  of  Attack 

Figure  21  shows  the  effect  of  an  elliptic  leading  edge  that  has  ten  degrees  of 
droop.  The  plot  represents  the  pressure  distribution  from  the  nose  of  the  leading  edge 
back,  and  again  suggests  that  the  flow  is  well  behaved  at  lower  angles  of  attack.  The 
bumps  in  the  pressure  distribution  after  the  pressure  spike  at  the  higher  angles  of  attack 
indicate  that  the  flow  may  be  separating.  However,  the  magnitude  of  these  pressure 
fluctuations  is  still  significantly  less  than  it  is  for  the  circular  leading  edge. 
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Figure  22.  Pressure  Distribution  (-Cp)  Over  Model  with  and  Ellipticity  Ratio  of  Four  and 

Various  Droop  Angles  at  Six  Degrees  Angle  of  Attack 

The  impacts  of  varying  droop  angle  can  be  seen  in  Figure  22  where  a  blade  model 
with  an  ellipticity  ratio  of  4  is  run  at  six  degrees  angle  of  attack.  For  the  case  where  there 
is  no  droop,  a  purely  elliptical  leading  edge,  the  pressure  spike  occurs  early,  and  has  a 
comparatively  high  value.  As  the  droop  angle  approaches  the  point  where  it  closely 
matches  the  incidence  angle,  the  magnitude  of  the  pressure  spike  is  the  lowest.  Further 
increases  in  the  droop  angle  past  the  angle  of  attack  result  in  the  value  of  the  pressure 
spike  rising.  Also,  as  the  droop  angle  is  increased,  the  point  along  the  blade  where  the 
pressure  spikes  occur  moves  back  away  from  the  leading  edge.  The  backside  of  the 
pressure  spike  shows  the  most  dramatic  decrease  in  the  case  where  the  droop  angle  is 
farthest  past  the  angle  of  attack,  suggesting  an  increase  in  the  magnitude  of  the  adverse 
pressure  gradient  and  the  formation  of  separation  bubbles. 
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2.  Boundary  Layer  Momentum  Thickness 

This  section  includes  results  for  a  selection  of  simulations  that  were  run.  The 
figures  presented  provide  a  generalization  of  the  behavior  of  the  boundary  layer.  The 
boundary  layer  momentum  thickness  was  measured  at  a  point  0.0350  meters  along  the 
body  of  the  blade,  approximately  forty  percent  of  the  chord  length.  At  this  point,  the  local 
Reynolds  number  was  517,619.  The  boundary  layer  momentum  thickness  was  arrived  at 
by  taking  velocity  data  in  the  boundary  layer  at  this  point,  and  was  reduced  according  to 
the  equation 


0  =  f  — 


'.-JL' 


V 


6/ 


dy 


(5) 


The  values  graphically  depicted  in  this  section,  and  those  tabulated  in  Appendix  B 
represent  the  sum  of  the  boundary  layer  momentum  thickness  on  the  pressure  side  (ps), 
and  suction  side  (ss)  of  the  model,  0SS  +0ps .  Figure  21  shows  an  example  of  the  boundary 


layers  on  the  suction  side  and  pressure  side  of  the  blade. 


Suction  Side  Pressure  Side 

Velocity  Profile  Velocity  Profile 


Figure  23.  Velocity  Profiles  on  the  Suction  Side  (ss)  and  Pressure  Side  (ps)  of  a  Plate  as 
Used  in  the  Calculation  of  the  Boundary  Layer  Momentum  Thickness,  6ss  +  0 
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Figure  24.  Circular  Leading-Edge  Boundary  Layer  Momentum  Thickness,  0SS  +  0  ,  at 

Various  Angles  of  Attack 

Returning  to  the  baseline  geometry  of  a  circular  leading-edge  shape,  Figure  24 
depicts  the  impact  of  angle  of  attack  on  the  boundary  layer  momentum  thickness.  In  the 
case  of  the  circular  leading  edge,  these  values  were  obtained  based  on  calculated  velocity 
profiles  taken  at  the  trailing  edge  of  the  model.  The  separation  that  occurred  on  the 
circular  leading-edge  model  made  calculating  the  boundary  layer  momentum  thickness  at 
the  original  calculation  point  unviable.  However,  for  the  undrooped,  zero  incidence 
simulations,  measurements  made  father  down  the  blade  where  the  flow  had  reattached 
allowed  for  this  calculation.  From  Figure  24,  it  can  be  seen  that  the  boundary  layer 
momentum  thickness  increases  as  angle  of  attack  increases.  Clearly  the  thickness  is  not 
linearly  dependent  on  angle  of  attack,  and  the  increase  becomes  more  dramatic  at  higher 
incidences. 
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Figure  25.  Momentum  Thickness,  0SS  +  0  ,  at  the  End  of  the  Plate  for  Various  Ellipticity 

Models 


Determining  the  effect  of  leading-edge  ellipticity  on  aerodynamic  losses  was  done 
by  comparing  the  boundary  layer  momentum  thicknesses  as  depicted  in  Figure  25.  As  the 
figure  shows,  the  momentum  thickness  of  the  circular  leading  edge  is  several  times  that 
of  the  elliptic  leading-edge  models.  With  no  droop,  the  blade  model  that  experienced  the 
least  losses  had  an  ellipticity  ratio  of  three.  A  very  shallow  loss  bucket  can  be  seen  in  the 
momentum  thicknesses  of  the  blades  with  neighboring  ellipticity  ratios. 
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Figure  26.  Boundary  Layer  Momentum  Thickness,  6SS  +  0ps ,  for  All  Ellipticities  at  Various 

Angles  of  Attack  and  Zero  Degrees  Droop 


Again,  as  is  seen  in  Figure  26,  the  circular  leading-edge  model  experienced 
significantly  higher  losses  that  did  the  elliptic  leading-edge  models.  This  figure  compares 
the  effect  of  angle  of  attack  on  the  boundary  layer  momentum  thickness,  and 
demonstrates  that  the  thickness  increases  with  increasing  angle  of  attack.  Figure  26  does 
not  have  any  thickness  values  at  an  incidence  of  13  degrees  because  all  models 
experienced  separation  at  the  point  of  boundary  layer  momentum  thickness  calculations 
at  this  incidence.  Additionally,  only  the  models  with  an  ellipticity  ratio  of  3.5  and  higher 
had  attached  flow  at  an  incidence  of  ten  degrees.  The  blade  model  with  an  ellipticity  ratio 
of  3.5  showed  the  least  dependency  on  angle  of  attack,  having  the  tightest  grouping  of 
momentum  thicknesses,  and  only  separating  at  13  degrees  incidence.  The  blade  model 
with  an  ellipticity  ratio  of  5.5  experienced  slightly  higher  losses  that  those  of  the  less 
elliptic  blades.  This  can  partially  be  attributed  to  the  increased  surface  distance  the  flow 
has  to  travel,  which  results  in  a  slightly  higher  boundary  layer  momentum  thickness. 
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Figure  27.  Boundary  Layer  Momentum  Thickness,  6SS  +  0ps ,  for  Model  with  Leading-Edge 
Ellipticity  Ratio  of  2.5  and  All  Droop  Angles  at  Various  Angles  of  Attack 

Analyzing  the  effects  of  droop  on  the  blade  model  with  an  ellipticity  ratio  of  2.5, 
Figure  27  shows  that  the  boundary  layer  momentum  thickness  increases  with  incidence. 
Separation  was  again  experienced  at  an  angle  of  attack  of  13  degrees  for  all  droops 
except  the  six  degree  case.  At  the  lower  incidences,  six  degrees  for  example,  the  lowest 
losses  are  observed  for  the  models  where  the  droop  angle  is  just  slightly  higher  than  the 
incidence  angle.  In  this  case,  drooping  the  leading  edge  to  six  degrees,  or  even  10 
degrees,  achieved  lower  losses  than  no  droop  or  drooping  it  to  13  degrees. 
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Figure  28.  Boundary  Layer  Momentum  Thickness,  9SS  +  0 ,  for  Model  with  Ellipticity  Ratio 

of  Four  and  All  Droop  Angles  at  Various  Angles  of  Attack 


A  similar  trend  to  that  depicted  in  Figure  27  is  seen  in  Figure  28.  In  general,  the 
boundary  layer  momentum  thickness  increases  with  increasing  angle  of  attack.  Flowever, 
drooping  to  where  the  8  angle  more  closely  match  the  incidence  angle,  a  results  in  lower 
losses  than  an  undrooped  leading  edge  at  incidence. 
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Boundary  Layer  Momentum  Thickness,  6SS  +  0ps ,  for  Model  With  Ellipticity 
Ratio  of  5.5  and  All  Droop  Angles  at  Various  Angles  of  Attack 


At  the  highest  ellipticity  ratio  tested,  5.5,  it  can  be  seen  that  drooping  the  leading 
edge  still  is  mostly  beneficial  when  compared  to  the  undrooped  case.  However,  as  Figure 
29  exemplifies,  at  very  low  incidences,  drooping  the  leading  edge  can  result  in  higher 
losses  that  it  would  at  higher  angles  of  attack.  In  this  case,  it  can  be  seen  that  the  greatest 
reduction  in  losses  occurs  for  the  models  that  have  a  droop  angle  closely  matching  the 
incidence,  indicating  the  impact  of  aligning  the  leading-edge  direction  with  the  angle  of 
attack  to  reduce  losses. 
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Figure  30.  Boundary  Layer  Momentum  Thickness,  0SS  +  0  ,  for  Models  of  All  Ellipticity 

Ratios  with  6.6  Degrees  Droop  at  Various  Angles  of  Attack 

In  Figure  30,  the  effects  of  ellipticity  with  constant  droop  are  examined  at  various 
angles  of  attack.  In  the  cases  where  the  flow  was  attached  at  the  point  of  boundary  layer 
momentum  calculation,  the  circular  leading  edge  has  losses  several  times  those  of  the 
elliptic  leading  edges.  Also,  from  Figure  30,  it  can  be  seen  that  the  lowest  losses  occur 
near  the  point  where  the  incidence  most  closely  matches  the  droop  angle. 


B.  DISCUSSION  OF  FLOW  FIELDS 

The  results  provided  give  insight  into  the  impact  of  the  leading-edge  geometry  on 
the  pressure  distribution  and  boundary  layer  momentum  thickness.  Analyzing  these 
results  allows  for  the  indication  of  boundary  layer  control  as  well  as  the  losses  that  are 
present  at  a  result  of  the  geometry  modification. 

Calculations  made  on  the  baseline  geometry  of  a  flat  plate  with  a  circular  leading 
edge,  allowed  for  the  determination  of  whether  or  not  a  particular  ellipticity  and  droop 
angle  were  advantageous  at  a  certain  angle  of  attack.  The  circular  leading-edge  model 
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experienced  heavy  separation  at  the  leading  edge,  which  resulted  in  poorly  developed 
flow  over  the  remaining  surface  of  the  blade  on  both  the  suction  and  pressure  sides. 

Even  where  the  flow  reattached,  a  separation  bubble  fonned.  Figures  31  and  32 
show  a  close  up  of  these  separation  bubbles  on  the  suction  side  and  pressure  side, 
respectively,  where  the  velocity  vectors  are  represented  by  arrows.  As  the  figures  show, 
in  the  separation  bubble  the  large  adverse  pressure  gradient  results  in  the  flow  actually 
traveling  backwards. 


Figure  3 1 .  Feading-Edge  Suction  Side  Separation  Bubble 


Figure  32.  Feading-Edge  Pressure  Side  Separation  Bubble 
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A  blade  with  a  higher  ellipticity  ratio  allowed  for  better  control  of  the  flow  over 
the  leading  edge,  as  seen  in  Figure  33.  Looking  at  the  flow  over  the  suction  side  of  the 
leading  edge  shows  that  even  at  an  angle  of  attack  of  ten  degrees,  the  flow  is  better 
controlled  than  in  the  case  of  the  circular  leading  edge.  However,  at  the  higher  angles  of 
attack,  such  as  shown  in  Figure  33,  the  flow  still  separates  and  transitions  to  turbulence, 
however  the  region  of  separated  flow  is  not  as  dramatic  as  for  the  circular  leading-edge 
model.  In  fact,  the  only  models  for  which  no  separation  was  noticeable  over  the  entire 
blade  body  had  ellipticity  ratios  of  four  and  5.5  with  zero  degrees  droop,  run  at  zero 
degrees  incidence. 


Figure  33.  Leading-Edge  View  of  Blade  Model  with  Ellipticity  Ratio  of  4,  Zero  Degrees 

Droop,  Ten  Degrees  Incidence 
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Figure  34.  Leading-Edge  View  of  Blade  Model  with  Ellipticity  Ratio  of  4,  Ten  Degrees 

Droop,  Zero  Degrees  Incidence 

Despite  the  apparent  advantages  of  the  drooped,  elliptic  leading  edge,  there  was  a 
more  distinct  impact  on  the  flow  over  the  pressure  side  of  the  airfoil  at  lower  angles  of 
attack.  Figure  34  shows  that  when  the  droop  angle  is  far  greater  than  the  incidence,  the 
velocity  boundary  layer  grows  much  quicker  on  the  pressure  side  than  it  did  for  the 
undrooped  case.  This  actually  results  in  a  separation  bubble  forming  at  the  point  where 
the  leading  edge  and  the  blade  body  are  blended  on  the  pressure  side.  This  separation 
bubble  decreases  in  size  and  eventually  vanishes  as  the  droop  angle  is  brought  closer  to 
the  angle  if  incidence.  On  the  suction  side,  the  separation  bubbles  also  decreased  in  size 
as  incidence  more  closely  match  the  droop  angle. 

The  results  presented  for  the  boundary  layer  momentum  thickness  are  more 
valuable  qualitatively  than  quantitatively.  The  prediction  of  transition  from  laminar  to 
turbulent  flow  still  remains  difficult,  and  the  tools  necessary  to  predict  this  transition  as 
applied  to  the  cases  run  are  not  available.  The  exact  values  arrived  at  for  boundary  layer 
momentum  thickness  hold  less  importance  than  the  trend  of  these  values.  The  main 
benefit  from  this  analysis  is  to  see  the  effects  of  leading-edge  geometry  modification.  A 
relative  comparison  of  the  momentum  thickness  between  the  different  cases  simulated 
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allows  for  the  impact  of  leading-edge  ellipticity  and  droop  to  be  observed.  Appendix  A 
provides  a  listing  of  all  cases  simulated,  and  whether  or  not  separation  occurred  on  the 
suction  side  and  pressure  side  near  the  leading  edge.  Also,  it  is  indicated  whether  or  not 
there  was  separation  at  the  point  of  velocity  measurements  used  for  the  calculation  of 
boundary  layer  momentum  thickness.  Appendix  B  provides  boundary  layer  momentum 
thickness  values  for  all  cases.  These  calculations  were  made  at  a  consistent  point  along 
the  plate,  however  for  comparison  purposes,  it  was  necessary  in  certain  instances  to 
provide  this  calculation  at  the  end  of  the  model,  or  the  trailing  edge,  as  indicated  by  a 
(TE). 


46 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  the  numerical  simulations  performed,  the  following  conclusions  were 

drawn: 

1 .  The  performance  and  ability  of  ESI-CFD  to  provide  results  acceptably 
close  to  experimental  values  was  verified.  The  pressure  distribution 
over  the  Walraevens  Cumpsty  model  provided  a  trend  that  was  as  close 
as  expected  to  the  experimental  results. 

2.  Numerical  results  for  all  simulations  were  obtained.  Pressure 
distributions  were  calculated  over  the  surface  of  each  model  that 
indicated  boundary  layer  behavior.  The  boundary  layer  momentum 
thickness  was  calculated  at  a  specified  point  on  each  model.  In  certain 
cases  the  point  where  the  boundary  layer  momentum  thickness  was 
calculated  had  to  be  changed  due  to  separation  occurring  at  the  initial 
calculation  point. 

3.  Increasing  the  incidence  that  the  models  were  simulated  at  generally 
resulted  in  a  greater  pressure  spike,  a  greater  adverse  pressure  gradient, 
and  an  increase  in  the  likelihood  of  separation  and  transition  to 
turbulence.  Increasing  the  ellipticity  allowed  for  smoother  acceleration 
and  deceleration  of  the  flow  over  the  leading  edge.  The  greater  the 
ellipticity,  the  less  likely  a  separation  bubble  would  occur.  Drooping  the 
leading  edge  to  match  the  angle  of  attack  resulted  in  the  greatest 
reduction  of  losses.  These  reductions  were  not  as  substantial  when  the 
droop  angle  was  less  than  or  greater  than  the  incidence  angle. 

4.  An  optimum  ellipticity  ratio  lies  somewhere  in  the  2.5  to  5.5  range.  The 
model  with  an  ellipticity  ratio  of  3  achieved  the  lowest  losses  at  zero 
degrees  incidence;  however  surrounding  ellipticities  achieved  boundary 
layer  momentum  thicknesses  that  were  very  close.  Over  a  range  of 
incidences,  the  model  with  an  ellipticity  ratio  of  3.5  displayed  the 
lowest  losses.  An  optimum  droop  angle  does  not  exist  over  a  range  of 
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angle  of  attack.  Rather,  an  optimum  droop  angle  exists  where  the 
leading  edge  is  drooped  to  most  closely  match  the  expected  incidence. 
A  moderate  droop  of  6.66  degrees  did,  however,  show  a  reduction  in 
losses  over  the  incidence  range  when  compared  to  a  circular  leading 
edge. 
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APPENDIX  A.  FLOW  SEPARATION 


Table  5  lists  all  the  regimes  tested,  and  are  marked  as  to  whether  or  not  the  flow 
separated  on  the  blade.  An  X  in  the  third  column  denotes  whether  or  not  there  was 
leading-edge  suction  side  separation.  The  fourth  column  denotes  whether  or  not  the  flow 
separated  on  the  pressure  side  of  the  leading  edge.  An  X  in  the  sixth  column  indicates 
that  the  flow  was  separated  where  velocity  measurements  were  made  on  the  suction  side 
of  the  blade  to  measure  the  boundary  layer  momentum  thickness.  This  indicates  that 
measurements  taken  at  this  point  were  not  used  in  the  calculation  of  the  momentum 
boundary  layer  thickness.  For  undrooped,  zero  incidence  cases  where  this  occurred,  the 
momentum  thickness  was  recalculated  at  the  trailing  edge  for  comparison  reasons. 


Table  5.  Simulation  Matrix  Indicating  Separated  Flows  Schemes 


Ellipticity 

Ratio 

Droop  Angle 

Incidence 

Suction  Side 

Leading-Edge 

Separation 

Pressure  Side 
Leading-Edge 
Separation 

Separated  Flow 
at  Point  of  9 
Measurement 

0.625 

0 

0 

X 

X 

X 

0.625 

3 

X 

X 

X 

0.625 

6 

X 

X 

X 

0.625 

10 

X 

X 

X 

0.625 

13 

X 

1 

0 

0 

X 

X 

X 

1 

3 

X 

X 

1 

6 

X 

X 

1 

10 

X 

X 

X 

1 

13 

X 

2.5 

0 

0 

X 

3 

X 

6 

X 

10 

X 

X 

13 

X 

X 

2.5 

3.333 

0 

X 

X 

3 

X 

X 

6 

X 

X 

10 

X 

X 

13 

X 

X 
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APPENDIX  B.  TABULATED  CALCULATIONS 


Table  6. 


Boundary  Layer  Momentum  Thicknesses 


Ellipticity 

Ratio 

Droop 

Angle 

(deg) 

Incidence 

(deg) 

Momentum  Thickness 
(meters) 

0.625 

0 

0 

0.00013992659014 (TE) 

0.625 

3 

0.625 

6 

0.625 

10 

0.625 

13 

0.00212813723579 

1 

0 

0 

1 

3 

0.00028797913689 

1 

6 

0.00036587474820 

1 

10 

1 

13 

1 

0 

0 

0.000460385996270825 (TE) 

1 

3 

0.00087850425430968 (TE) 

1 

6 

0.00103134429584577 (TE) 

1 

10 

0.00157126980558499 (TE) 

1 

13 

0.00198567433925469 (TE) 

2.5 

0 

0 

0.000100094059520359 

3 

0.000103623994917905 

6 

0.000154502573900588 

10 

13 

2.5 

0 

0 

0.000140104356256184 (TE) 

3 

0.000296696428376302 (TE) 

6 

0.000539087536161559 (TE) 

10 

0.00111492335987182 (TE) 

13 

0.0017212911378331 (TE) 

2.5 

3.333 

0 

0.000111940154827703 

3 

0.000107007881792047 

6 

0.000132018572754849 

10 

0.000393112974704518 

13 

2.5 

6.666 

0 

0.00010999359999889 

3 

0.000106563326027226 

6 

0.000131588168817812 

10 

0.000361654408114119 
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13 

0.000580175862682864 

2.5 

10 

0 

0.000135036565426467 

3 

0.000112853954569468 

6 

0.000128443059324597 

10 

0.000408886437067048 

13 

2.5 

13.333 

0 

0.000170979688948428 

3 

0.000146508049516922 

6 

0.000157988001399255 

10 

13 

3 

0 

0 

9.62451 589779253e-005 

3 

0.000103487255926937 

6 

0.000123698773002648 

10 

13 

3 

0 

0 

0.000137071890747477 (TE) 

3 

0.000143351992026842 (TE) 

6 

0.000312874597826619 (TE) 

10 

0.000609292125163118 (TE) 

13 

0.00114217567150724 (TE) 

3.5 

0 

0 

9.6466541 6381 034e-005 
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APPENDIX  C.  PRESSURE  DISTRIBUTION  PLOTS 


Figure  35.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  2.5  and  Droop  Angle  of  3.33  degrees 
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Figure  36.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  2.5  and  Droop  Angle  of  6.66  degrees 


Figure  37.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  2.5  and  Droop  Angle  of  10.0  degrees 
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Figure  38.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  2.5  and  Droop  Angle  of  13.33  degrees 


Figure  39.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 

Ellipticity  Ratio  of  4.0  and  Droop  Angle  of  0.0  degrees 
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Figure  40.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  4.0  and  Droop  Angle  of  3.33  degrees 


Figure  41.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  4.0  and  Droop  Angle  of  6.66  degrees 
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Figure  42.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  4.0  and  Droop  Angle  of  10.0  degrees 


Figure  43.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  4.0  and  Droop  Angle  of  13.33  degrees 
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Figure  44.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 

Ellipticity  Ratio  of  5.5  and  Droop  Angle  of  0.0  degrees 


Figure  45.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  5.5  and  Droop  Angle  of  3.33  degrees 
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Figure  46.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  5.5  and  Droop  Angle  of  6.66  degrees 


Figure  47.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  5.5  and  Droop  Angle  of  10.0  degrees 
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Figure  48.  Pressure  Coefficient  (-Cp)  over  Model  from  Geometric  Leading  Edge  with 
Ellipticity  Ratio  of  5.5  and  Droop  Angle  of  13.3  degrees 
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APPENDIX  D.  MATLAB  CODE 


A.  LEADING-EDGE  BLENDING  CODE 

The  following  code  was  used  to  calculate  the  leading  edge  and  plate  geometries 
used  in  the  simulations. 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%Jonathan  D.  Powell 
%Thesis  Work 

%This  File  Started  6  Mar  2005 

%Flat  Plate  with  a  Modified  Leading  Edge 

%Establish  the  flat  plat 

%  We'll  make  this  flat  plate  approximately  8%  thick  (excluding  leading 

%  edge) 

clc 

fp_x=0:0. 5:100; 
for  i=l  :length(fp_x) 
fp_suct_y(i)=4; 
fp_press_y(i)=-4; 
end 

%plot(fp_x,fp_suct_y,fp_x,fp_press_y) 

%axis([-20  110-30  30]) 

%axes  equal 

%Now  to  make  modifications  to  the  leading  edge 
%Starting  with  varying  eccetricity  of  unrotated  ellipse 
b=fp_suct_y(l); 

%incr=l:0.5:5.5; 
incr=[l  2.5  4  5.5]; 
incr=incr(2) 
a=b.*incr; 
start=length(incr)/4; 
stop=(3/4)  *  length(incr) ; 
for  i=l  :length(incr) 

%for  j=l  :length(incr) 
ecc(i)=sqrt(l-((b)A2)./((a(i)).A2)); 

[coord_x,coord_y]  =  ellipse  1  (0,0, [a(i)  ecc(i)],0); 
end 

qcc=qcck, 
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%rot=0: 1 5/(length(incr)- 1 ):  1 5; 
rot=[0  3.33333  6.66667  10.0  13.33333]; 
rot=rot(5) 
for  j=l  :length(incr) 
figure 

for  k=l :  length!  rot) 

[coordxrot.coordyrot]  =  ellipsel(5,0,[a(j)  ecc(j)],rot(k),  [(75)  (330)]); 
%figure 

%subplot(3,2,k) 
for  tt=l:50 

if  ((coord_y_rot(tt)'  <  4)) 
y_tip_a(tt)=coord_y_rot(tt)'; 
else 

y_tip_a(tf)=fp_suct_y(l ); 
end 
end 

forpp=l:50 

if  (min(coord_y_rot)  <  -4)  &  (coord_y_rot(  1 00)'  >  -4) 
y_tip_b(pp)=coord_y_rot(pp+50)'; 
else 

y_tip_b(pp)=fp_press_y(  1 ); 
end 

y_tip=[y_tip_a  y_tip_b]; 
end 

forss=l:15 

if  y_tip_b(ss+35)  >  -4 
y_tip_b(  ss+3  5  )=fp_press_y(  1 ) ; 
else 

y_tip_b(ss+35)=coord_y_rot(ss+85); 

end 

end 

y_tip=[y_tip_a  ytipb]; 
for  xx=l:50 

if  fp  x(xx)  <  coord_x_rot(  1 ) 
fp_x_suct(xx)=coord_x_rot(  1 ); 
else 

fp_x_suct(xx)=fp_x(  xx) ; 
end 
end 

forrr=l:50 

if  fp_x(rr)<coord_x_rot(  100) 

fp_x_press(rr)=coord_x_rot(  100); 


else 
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fp_x_press(rr)=fp_x(rr); 

end 

end 

fp_x_suct_a=[fp_x_suct  fp_x(50:200)]; 
fp_x_press_a=[fp_x_press  fp_x(50:200)]; 

%  plot(fp_x,fp_suct_y,fp_x,fp_press_y,coord_x_rot',coord_y_rot,'-') 

%plot(fp  x  sucl  a. Ip  suet  y,fp  x  press  a.fp  press  _y, coord  x  rot'.y  tip.'.', coord  x  rct(l),y  tip(l),'*', coord  x  rod  100)', 

y_tip(100),'*',coord_x_rot,coord_y_rot,5,0,'*',x,y,x_new,y_new,'.') 

%axis  equal 
%axis([-10  20-10  20]) 

%axis  on 
end 
end 

x_lines=2:10; 
a=coord_x_rot(3 :5); 
b=y_tip(3:5)'; 
pa  =  polyfit(a,b,  1 ); 
m_a=p_a(  1 ); 
b=p_a(2); 

c=coord_x_rot(  14:16); 
d=y  tip(  14: 16)'; 
pc  =  polyfit(c,d,  1 ); 
m_c=p_c(l); 
b_c=p_c(2); 

yc=-l/m_c*(x_lines-coord_x_rot(15))+y_tip(  15); 
for  1=1:9 

x_c_sta(i)=coord_x_rot(4); 

end 

pp_c=polyfit(x_lines,y_c,  1 ); 

y_cc=pp_c(  1  )*(coord_x_rot(4))+pp_c(2); 

r_c=sqrt(  (coord_x_rot(  1 5  )-x_c_sta(  l))A2+(y_tip(  1 5  )-y_cc)A2); 

x_circ_c=coord_x_rot(  1 5):0.3  :coord_x_rot(4); 

ycircc=sqrt(r_cA2-(x_circ_c-x_c_sta(l)).A2)+y_cc; 

e=coord_x_rot(90:92); 

f=y_tip(90:92)'; 

p_e=polyfit(  e,f,  1 ); 

m_e=p_e(  1 ); 

b_e=p_e(2); 

g=coord_x_rot(79:8 1 ); 

g=coord_x_rot(79 :8 1 ); 

h=y _ tip(79:8 1 )'; 

p_g=polyfit(g,h,l); 

m_g=P_g(l); 
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b_g=P_g(2); 

y_g=-l/m_g*(x_lines-coord_x_rot(80))+y_tip(80); 
for  i=l:9 

x_g_sta(i)=coord_x_rot(9 1 ); 
end 

pp_g=polyfit(x_lines,y_g,  1 ); 
y_gg=pp_g(l  )*(coord_x_rot(9 1  ))+pp_g(2); 

%r_g=sqrt((coord_x_rot(78)-x_g_sta(l))A2+(y_tip(78)-y_gg)A2); 
r_g=sqrt((coord_x_rot(80)-x_g_sta(l))A2+(y_tip(80)-y_gg)A2); 
xcirc_g=coord_x_rot(80):0.3:coord_x_rot(91); 
y_circ_g=sqrt(r_gA2-(x_circ_g-x_g_sta(  1  )).A2)+y_gg; 
x  circ  LE=0:0. 1 :5 

y_circ_LE=sqrt(4A2-(x_circ_LE-5).A2) 

coord_x_rot=[coord_x_rot(l:3)'  fliplr(xcircc)  coord_x_rot(  16:77)'  x  circ  g  coord_x_rot(93:100)']; 
y  tip=[y  tip(  1:3)  fliplr(ycircc)  y  tip(  16:77)  y_circ_gy_tip(93:100)]; 

plot(fp_x_suct_a,fp_suct_y,fp_x_press_a,fp_press_y,coord_x_rot',y_tip,'.',coord_x_rot(l),y_tip(l),’*',coord_x_rot(length(c 
oord  x  rot))',y  tip(length(y  tip)),'*',x  circ  g.y  circ  g,'*',x  circ  c,y  circ  c,'*',x  lines.y  g,x  g  sta,y  g)%,x  lines, y  c,x  c  sta,y  c.x 
g_sta(l),y_gg,'*') 

%axis  equal 

%axis([-8  12-10  10]) 

axis([-10  25  -25  10]) 

%end 

%end 

%runs_matrix_LE=[coord_x_rot,  y  tip'] 
runs_matrix  fp_x=[fp_x_suct_a  fp_x_press_a]'; 
rans_matrix_fp_y=[fp_suct_y  fp_press_y]'; 
disp('X') 

aaa=[coord_x_rot  fp  x  suct  a  fp_x_press_a]' 
disp('Y') 

aab=[y_tip  fp  suct  y  fp_press_y]' 

%disp('FP  X') 

%aac=runs_matrix_fp_x 
%disp(  'FP  Y') 

%aad=runs_matrix_fp_y 

%runs_matrix_fp=[runs_matrix_fp_x  runs  matrix  I'p  y] 


B.  ELLIPTIC  CALCULATIONS  CODE 

The  MATLAB  code  below  was  used  to  calculate  the  ellipticity  ratios  of  the 
undrooped  models,  and  the  arc  length  of  their  leading  edges.  These  calculations  were 
made  based  on  the  distance  from  the  geometric  leading  edge,  at  the  semimajor  apex,  to 
the  point  where  the  elliptic  leading  edge  joined  the  plate  body,  at  the  semiminor  apex. 
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ER=[0.625  1  2.5  3  3.5  4  5.5] 

%Semiminor 

a=[0.004]; 

%  Semimajor 

b=[0.625*a  l*a  2.5*a  3*a  3.5*a  4*a  5.5*a]; 

R_min=(b.A2)./a 

R_maj=(a.A2)./b 

t=0.008 

s=(pi.*  [3  *(a+b)-sqrt((3  *a+b).  *(a+3  *b))])/4 
two_s_over_t=2  *s./t 

WalCump_circ_a=[0.0063]; 

WalCump_circ_b=[0.0063]; 
t_WalCump_circ=0.0 1 26 
WalCump_ellip_a=[0.0066]; 

WalCump_ellip_b=[0.0125]; 
t_WalCump_ellip=0.0 1 32 

s_WalCump_circ=(pi.  *  [3  *(WalCump_circ_a+WalCump_circ_b)- 

sqrt((3*WalCump_circ_a+WalCump_circ_b).*(WalCump_circ_a+3*WalCump_circ_b))])/4 

s_WalCump_ellip=(pi.*[3*(WalCump_ellip_a+WalCump_ellip_b)- 

sqrt((3*WalCump_ellip_a+WalCump_circ_b).*(WalCump_ellip_a+3*WalCump_ellip_b))])/4 
twosovertW  alCump_circ=2  *  sW  alCumpcirc  ./t_  W  alCumpcirc 
%for  elliptical  model,  in  order  to  get  a  2*s/t  of  20,  we  need  to  go  an 
%additional  18.4292*t/2 

extra_WalCump_circ=18.4292*t_WalCump_circ/2 
%blendpoint  for  the  circle  is  at  2s/t=l .5708 

two_s_over_t_WalCump_ellip=2  *  s_WalCump_ellip./t_WalCump_ellip 
%for  elliptical  model,  in  order  to  get  a  2*s/t  of  20,  we  need  to  go  an 
%additional  1 7.28 1 6*t/2 

extra_WalCump_ellip=l  7.28 1 6*t_WalCump_ellip/2 
%blendpoint  for  the  ellips  is  at  2s/t=2.7184 

plot(ER,s, xlabel('Arclength  on  Ellipse  from  Semimajor  Apex  to  Semiminor  Apex,  meters') 

ylabel('Ellipticity  Ratio') 

vec=[0  0  0  0  0  0  0] 

plotvec_x_l=[0,s(  1 )] 

plotvec_x_2=[0,s(2)] 

plotvec_x_3=[0,s(3)] 

plotvec_x_4=[0,s(4)] 

plotvec_x_5=[0,s(5)] 

plotvec_x_6=[0,s(6)] 

plotvec_x_7=[0,s(7)] 

plotvec_y_l=[R_maj  ( 1  ),R_min(  1 )] 

plotvec_y_2=[R_maj  (2),R_min(2)] 
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plotvec_y_3=[R_maj(3),R_min(3)] 
plotvec_y_4= [  Rmaj  (4),R_min(4)] 
plotvec_y_5=[R_maj(5),R_min(5)] 
plotvec_y_6=[R_maj  (6),R_min(6)] 
plotvec_y_7=[R_maj  (7),R_min(7)] 

plot(plotvec_x_  1  ,plotvec_y_  1  ,plotvec_x_2,plotvec_y_2,,-A',plotvec_x_3  ,plotvec_y_3 :',plotvec_x_4,plotvec_y_4,,— 
',plotvec_x_5,plotvec_y_5,'-.',plotvec_x_6,plotvec_y_6,,-*',plotvec_x_7,plotvec_y_7,'-o') 

%plot(vec,R_maj,,*,,s,R_min,'o') 

legend('ER=0.6257ER=r,'ER=2.5','ER=3','ER=3.57ER=47ER=5.5') 

xlabel('Arclength  From  Geometric  Leading  Edge  To  Blend  Point,  meters') 
ylabel('Radius  of  Curvature,  meters') 
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APPENDIX  E.  SUPPLEMENTAL  FIGURES 


Figure  49  represents  the  variation  of  leading-edge  arclength  with  ellipticity  ratio. 
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Figure  49.  Variation  of  Elliptical  Arclength  with  Ellipticity  Ratio 
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